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ABSTRACT 

We describe the latest release of AtomDB, version 2.0.2, a database of atomic data and a plasma 
modeling code with a focus on X-ray astronomy. This release includes several major updates to the 
fundamental atomic structure and process data held within AtomDB, incorporating new ionization 
balance data, state-selective recombination data, and updated coUisional excitation data for many 
ions, including the iron L-shell ions from Fe +16 to Fe +23 and all of the hydrogen- and helium-like 
sequences. We also describe some of the effects that these changes have on calculated emission and 
diagnostic line ratios, such as changes in the temperature implied by the He-like G-ratios of up to a 
factor of 2. 



1. INTRODUCTION 

X-ray spectra from astrophysical sources reveal both 
the constituent elements of those objects and the physics 
occurring within them. Successfully extracting this infor- 
mation from spectra requires both a model of the plasma 
and a large collection of data detailing the various atomic 
processes occurring within that plasma. 

The recognition of the importance of the dielec- 
tronic recombinatio n (DR) process for plasmas, even at 
high temperatures (Burgcsa 119641 ) . led to huge strides 
in modeling collisionally ionized, optica l ly-thi n astro- 
physical plasmas. The ICox fe Tucker! (|1969l) collec- 
tion was one of the earliest attempts to collect atomic 
data for these models. Since then there have been 
a steady stream of refinements to both the coUisional 
plasma models and the atomic data which underpin 
them. Differe nt evolutions of s uch models have in - 
cluded those of ICox fc Daltabuitl Jl97lh. iMewel (fl972T) 
lLandini fc Monsignori Fossil ( 19721) . iRavmond fc Smith! 
(|1977t ). and lBrickhouse et all (|199a ). 

Both the quantity and quality of relevant atomic data 
continue to grow: increases in computational power avail- 
able allow for improved calculations, while advances in 
experimental methods and equipment allow more ac- 
curate measurements of a wider range of quantities. 
In modern analysis of collisionally ionized X-ray as- 
trophysical spectra the re are three widely- used atomic 
dat abases: SPEX v2 (jKaastra et al.l ll996l). CHIANTI 
v7 (lLandi et al.ll2012[ ) and AtomDB vl.3.1 jSmith et al.l 
l2001al) . Each database has a slightly different focus: 
CHIANTI's main focus is on the EUV wavelengths for 
analyzing solar spectra, while SPEX and AtomDB focus 
on the X-ray ranges. Each of these databases continues 
to undergo periodic review and is updated as newer data 
become available; this paper introduces the release of 
AtomDB v2.0.2, describing the new data and improve- 
ments, how they have been implemented and what is 
planned for future releases. 
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2. CALCULATION AND SELECTION OF ATOMIC 

DATA 

2.1. Priorities and Techniques 

In this paper we primarily discuss the theoretical 
sources of atomic data we have used for AtomDB v2.0.2. 
We will therefore begin by outlining a few of the guid- 
ing principles in our collection of data, and then briefly 
outlining the salient points of each major theoretical 
method. This is not intended as a comprehensive re- 
view of the atom ic physics; for such a do cument see a 
text book such as iBovle fc Pindzolal (|2005f ). 

2.1.1. Theoretical vs Experimental Data 

Accurate calculation of the emission from a 
collisionally-ionized plasma requires a broad range 
of atomic data, both in terms of the number of ions 
and processes included. The fundamental atomic data 
required can be split into structural and coUisional 
categories. Structure data include the atomic energy 
levels and the radiative transition rates among them. 
CoUisional data describe interactions between the ion 
and another particle in the plasma - most commonly, 
an electron, although proton impact can have an effect 
on low-lying levels. The simplest atomic processes that 
can arise from collisions are electron impact excitation 
(EIE), electron impact ionization (EII), and radiative 
recombination (RR) . More complex multi-step processes 
also occur, such as excitation-autoionization (EA) and 
dielectronic recombination (DR). 

The ideal method for obtaining an accurate atomic 
database would be to have a large and accurate collec- 
tion of experimental data covering all of these processes. 
In some cases, particularly for transiti on wavelengths, 
this c an be done: the NIST database (jRalchenko et al.l 
1201 1|) collects wavelengths of line emission from a wide 
range of ions, from which highly accurate energy lev- 
els are calculated along with uncertainties. However for 
many processes, experimental measurements with suffi- 
cient breadth of coverage are unavailable: measuring col- 
lision strengths for one ion can take many months, while 
investigating it over the wide range of energies of astro- 
physical interest can take longer, or require a different 
apparatus. In addition, systematic errors can affect ex- 
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perimental measurement and are sometimes very difficult 
to remove. Finally, sometimes identifying what has been 
measured is not simple: a line with a wavelength may 
be observed, but attributing that to the correct energy 
levels requires some theoretical work and can be difficult 
in very complex ions. Therefore it is customary to use 
theoretical calculations in databases, with the available 
experimental information employed as a quality check. 

One exception to this is line wavelengths: these can be 
measured with a high degree of accuracy. When using a 
fitting routine to match observed and modeled spectra, 
a small shift in wavelength can badly affect the results, 
inflating the fit statistic and occasionally causing other 
ion abundances to be affected. For this reason, observed 
wavelengths are preferred to the calculated wavelengths 
where they are available and the upper and lower levels 
can be accurately identified. 

2.1.2. Plasma Timescales 

The various processes in plasma occur on different 
timescale s. A full discus s ion of these timescales is pre- 
sented in lSummers et all ()2006t ). but it is sufficient here 
to note that r g ~ r TO 3> r D ^> r a 3> r ee , where timescales 
for change are: r g - the ground state ion population (typ- 
ically due to ionization or recombination processes); r m 

- metastable states (those without a fast radiative decay 
to a lower level) ; t - ordinary excited states (those with 
energy lower than the ionization potential of the ion); r a - 
auto-ionizing states (those with energy above the ioniza- 
tion potential, which can ionize spontaneously); and T ee 

- the thermalization timescale for electron-electron colli- 
sions to create a Maxwellian energy distribution in the 
plasma. A common assumption for astrophysical plas- 
mas, and one that is made in the AtomDB database, is 
that the electrons have a Maxwellian energy distribution: 
Maxwellian-averaged collision strengths are stored for all 
ions. 

2.1.3. Electron Impact Excitation 

Much of the new data compiled in v2.0.2 contains im- 
provements in the EIE calculation. A variety of methods 
exist for calculating collisional excitation by electrons, 
in particular the distorted wave (DW) and i?-matrix ap- 
proaches. The DW technique is a perturbative approach, 
in which the electron is modeled as a wave perturbing 
the N or N+l electron system. In general, this approach 
works best for highly ionized targets (more than 2 or 3 
times ionized), where the nucleus-electron potential can 
be treated as the dominant term, with all other interac- 
tions as lesser perturbations. It includes initial and final 
states of the system, but omits the many intermediate 
resonance states. These can be subsequently added back 
in in an isolated resonance approximation (i.e. assuming 
there is no correlation between each resonance). 

Close-coupled techniques, of which i?-matrix 
(|Burke et al.lll97ll ) is the most widely used for collisional 
excitation modeling, handle resonances naturally. In 
this model, the calculation is split into three regions 
based on the incoming electron-nuclear distance. The 
most significant change is in the inner region: the N+l 
electron system wave functions are constructed from N 
and N+l electron channels, with coupling possible with 
all the possible final states for the scattered electron. 
This allows modeling of resonances not as a separate 



process, but as an integral part of the calculation. It 
also allows i?-matrix to be used to consistently provide 
ionization and, to some extent, recombination data. 
The inclusion of resonances, and some damping terms 
required by the technique to avoid over counting, enables 
the i?-matrix method to give a more accurate treatment 
of lower energy collisions. It is, however, also a much 
more computationally intensive technique. For example, 
i?-matrix calculations of H- and He-like systems (those 
with 1 and 2 electrons, respectively) rarely include levels 
above the n — 6 shell, while DW calculations can include 
much higher levels with ease. There are some complex 
ions which are simply not possible to model well with 
i?-matrix, though these ions are generally heavier than 
those of astrophysical interest. In the end, a trade-off 
is required between the accuracy of i?-matrix and the 
wide coverage of ions and transitions offered by DW. 

The i?-matrix package has been modified and extended 
to include a variety of different models since its origi- 
nal introduction. The energy levels of heavier elements 
must be split in to their J-resolved sub levels for accu- 
rate spectroscopy. The Intermediate Coor dinate Frame 
Transformation (ICFT, iGriffin et a0ll998H method pre- 
forms much of the calculation in LS coupling, then ap- 
plies a transform to the results to obtain IC coupled (and 
therefore J-resolved) data. This is beneficial due to the 
much lower number of energy terms in an LS coupled 
calculation compared to a similar IC problem; the com- 
putation time scales as approximately n 3 , where n is the 
number of levels or terms involved. 

In Breit-Pauli i?-matrix, outlined in iHummer et al.l 
( 1993), the Hamiltonian for the system is formed by sim- 
plifying the fully-relativistic Dirac equation by treating 
the relativistic corrections as perturbations. For heavier 
elements (Z > 36), the fully relativistic Dirac equ ation 
should be solved: the DARC (jAit-Tahar et al.lll996f ) vari- 
ant of i?-matrix handles this, though it is usually not nec- 
essary for transitions from astrophysically relevant ions 
in the X-ray bandpass. Data from both the ICFT and 
Breit-Pauli methods have been included in this release. 

For ions with low ionization energies, in particular 
near- neutral ions, loss of excited electrons to the con- 
tinuum (i.e. ionization) can be a significant process 
and difficult to separate from excitation processes with 
near-ionization threshold energies. The introduction of 
pseudo-states to represent the continuum states signifi- 
cantly improves collision cross section calculations, usu- 
ally by reducing excessively high values. The i?-matrix 
with Pseudo States (RMPS) code (jBartschat et al.lll996f) 
incorporates this, along with its relativistic counterpart 
DRMPS. 

2.2. Selection of Atomic Data 

In selecting the data for inclusion in this version of 
AtomDB, we have focused most attention on those areas 
where there was a known problem. For example if line 
ratios were consistently not matching observations, and 
no other explanation was satisfactory, we attempted to 
obtain updated rates. Hence we have focused on the He- 
like and Fe L-shell ions, where discrepancies have been 
noticed, and in extending H- and He-like ions to higher 
n-shells. The goal of the AtomDB project is not to gen- 
erate atomic data, but to collect it from existing calcu- 
lations and distribute the data to the X-ray astronomy 
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community in useful formats. When there is no readily 
available published atomic data for a particular process 
or transition we will calculate and include it, however we 
will replace this calculated data with published results 
later when possible. 

Finding multiple calculations of the same rates is rel- 
atively rare, but it docs require that a decision be made 
as to which to include. In many cases, there is a clear 
improvement in the technique used in one of the calcu- 
lations which means the results are likely to be signif- 
icantly better than other data sets, such as comparing 
an i?-matrix collision calculation with a DW one with 
obvious effects on the lower energy collision strengths. 
However, it is not as simple as using the latest data - us- 
ing a newer technique for a calculation does not matter 
if it is not done well. We therefore look for comparisons 
with observation or laboratory measurements to select 
the best technique. Often multiple calculations have sim- 
ilar results, which makes the choice less relevant. There 
are relatively few cases where the data from two sources 
diverge strongly and there is little experimental data to 
provide any guidance; when these circumstances arise, 
we are forced to make a choice. We then follow the liter- 
ature post-release to discover if there are any recurring 
issues. 

3. NEW DATA 

Changes have been made throughout AtomDB. Several 
ions which were not previously included are now included 
in the database, while data for nearly all other ions have 
been improved by including more accurate calculations, 
by increasing the number of energy levels and transi- 
tions included for the ion, or a combination of both. In 
addition, major changes have been made to the ioniza- 
tion and, particularly, recombination data, with some 
significant changes to the resulting ionization balance 
and therefore to line emissivities. Here we outline each 
change in detail; in section [4] we will summarize the ef- 
fect of these changes on the spectra observed by users. 
We note that in this release there are no changes to the 
format of the data fil es in which the datab ase is stored; 
these are described in ISmith et al.l (|2001bfl . 

3.1. Ionization and Recombination Rates 

As experimental and theoretical methods improve, 
there are occasional efforts to collect the rates for ion- 
ization and recombination processes, including EII, EA, 
RR and DR recombination, for easy use by astrophysi- 
cists. Previous examples in clude the widely-used r ates o f 
lArnaud fc Raymond! (|1992l ) and IMazzotta et al.l (119981): 
the la t ter wa s included in AtomDB vl.3.1. IBrvans et al.l 
(|2006l I20Q9D have produced a new compilation for ele- 
ments from hydrogen [Z = 1) to zinc (Z=30), which we 
have included in AtomDB v2.0.2. This provides a com- 
plete set of new rates for ionization and recombination 
for all astrophysically relevant ions from fully stripped up 
to Na-like systems (where Na-likc refers to the recombin- 
ing ion). Both the rates and an equilibrium ionization 
balance table are included in AtomDB v2.0 .2. 

In their papers IBrvans et all (|2006l [2009) compare the 
electron-impact ionization (EII) and recombination data 
with earlier data sets and with experimental measure- 
ments. The most significant changes, particularly for 



TABLE 1 

Data sources for the DR 

STATE-SELECTIVE RATES. 



Sequence 



Reference 



H-like Badncll (2006b) 

He-like Bautista & Badncll (2007 

Li-like Colgan et al. (2004) 

Be-like Colgan et al. (2003) 

B-like Altun et al~72004T 

C-like Zatsarinnv et al. (2004b') 

N-like Mitnik fc Badnell (20041 

O-like Zatsarinnv et al. (2003) 

F-like Zatsarinnv et al. (2006) 

Ne-like Zatsarinnv et al. Q004a) 

Na-like Altun et al. (2006) 

Mg-like Altun et al. (2007) 



many L-she ll ions, have been in t he DR rates from the 
widely used lMazzotta et all ( 1998^) data. Cal culations of 
PR r ate coefficients bylB^dnelj ([2006a|) an d iGul (f2 003a. 
I2004D agree to within 35% of each other. IBrvans et al.l 
(2006) also find that the DR rates agree to within 35% 
for strong resonances when compared with results from 
storage ring experiments (see their paper for references) . 
The exception to this is An = core transitions at low 
temperature (T < 35,000-fQ, where results diverge: in 
this temperature range single and doubly ionized ions are 
dominant and their populations should be treated with 
care, while at higher temperature no such complications 
arise. Due to the gene ral good agreem ent of the meth- 
ods, the PR data fr om iBadnelll (|2006ah are used in the 
IBrvans et al.l (|2006| ) compilation where possible as they 
cover a wider range of ions (the H-like through Na-like 
systems fr om He through to Zn). The remainder is from 
the earlier Mazzott a et al.l (119981) data. The data s et for 
PR into nitrogen-like ions ([Mitnik fc Badnelf 2004) have 
been revised by their authors since their publication and 
we have included this revised data here. The data sources 
for each isoelectronic sequence are listed in Tab le [J 

RR r ates have been taken from calculations by I Badnell 
(2006 a[). T hese have been compa r ed wi th values from 
IGul (|2003H ) and IVerner fc Ferlandl ([1996), and agree to 
within 10% with the former, and 5 % w i th the latter. 
There is a larger divergence with the |G3 (I2003bl ) data, 
up to 20%, for H-like ions at temperatures in the upper 
range of their collisional ionization window, but in these 
cases the total recombination rate is dominated by the 
PR process so this effect is not important. Again due to 
the broader range of ions covered (bare nucleus to Na- 
like i ons from H to Zn) the Badnell data set is adopted, 
with IMazzotta et al.l [1998) data for the remaining ions. 

For EII the re are wider disagreements among data col- 
lections. The lPerel (]2007l ) study collected ionization data 
from all ions of all elements from H to Zn, from both lab- 
oratory and theoret ical sources. Another collection by 
iMattioli et al.l (|2007t ) covered a large subset of these ions: 
significant disagreement was fo und for many ions, with 
up to a factor o f four difference. IBrvans et al.l ([20091 ) use 
the iPerel (|2007f ) data due to its broader scope, but the 
authors note (and we agree) that this discrepancy should 
be revisited and if possible resolved in the future. 

3.2. Final State Resolved DR & RR Rates 
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The DR and RR data in the iBrvans et all (|2006l ) col- 
lection contain total rate coefficients for recombination 
from the ground state of one ion to the next. The data 
on which these are based (Table [1]) are resolved into re- 
combination into each excited state of the recombined 
ion. Inclusion of these rates when calculating the excited- 
state populations of the recombined ion will affect emis- 
sion lines and diagnostic line ratios; they will do so even 
more in plasmas dominated by recombination. For this 
reason we have included state-selective data for recom- 
bination to all H-like, He-like and Li- like ions, and all 
iron ion s from Fe +16 to Fe +22 . As provided by IBadnell 
(20063), DR & RR source data are split into level- 
resolved rates for capture into lower n (typically n < 8), 
with capture rates into higher n provided as totals into 
each n-shell. Since the levels identified in the DR & RR 
calculations and those in AtomDB are rarely identical, 
some level matching and cascade calculations have been 
required. 

Using the Atjtostructure (]Badnelllll986r ) code, tran- 
sition probabilities have been calculated for all the en- 
ergy levels which appear in either the AtomDB data 
sets or the DR/RR calculations. In addition, transition 
probabilities between n shells have been estimated using 
the h ydrogenic approximation of Burg ess fe Summers! 
(1976). Using this data, projection matrices have been 
formed to project the capture into the high-n levels onto 
the lower-n levels which are included in AtomDB. During 
the cascade, the electron density is assumed to be neg- 
ligible and thus there is no collisional re-ionization/re- 
distribution. This cascaded population is then dis- 
tributed according to the statistical weight of each level 
within the n shell. The direct rate into each of the lower 
levels is then added to cascaded contribution, to obtain 
a total effective recombination rate to each level in the 
AtomDB database. 

This treatment of the cascade process is not ideal: 
firstly, the photons emitted during the cascade are not 
tracked and therefore are not included in any spectral 
models or cooling function estimates. For the purposes 
of AtomDB, this is relatively unimportant since most cas- 
cade emission will not be in the X-ray wavelengths, and 
the total omitted radiated power is small. In addition, 
the treatment of DR satellite lines is not consistent with 
these data: in AtomDB satellite lines are tabulated as 
a function of temperature obtained from separate calcu- 
lations. The APEC code tracks the total DR rate from 
each approach to look for large discrepancies, with the 
average being a factor of 2 for strong transitions, but 
there are outliers which are much larger. They would 
be better handled in a self-consistent manner during the 
original calculation of the DR rates. 

3.3. H-like Ions 

The data for hydrogenic ions have been upgraded in 
two major ways: (1) they have been extended to include 
higher n-shells (from n = 5 to n = 10) and (2) where 
available, data of a higher quality have been used for 
ions of astrophysical interest. 

The existing data in AtomDB vl.3.1 for H-like 
ions come from a combination of the fully re l ativis - 
tic Dirac R- matrix calculations of iKisielius et al.l (|1996l ) 
for excitation from the ground, first excited and sec- 
ond excited state, and the non-relativistic calcula- 



TABLE 2 

The sources of data and maximum n shells for 

i?-MATRIX COLLISIONAL EXCITATION CALCULATIONS 
OF H-LIKE IONS. 



Ion 


Reference 




Coupling 


H° 


Anderson et al. 


(2002) 


5 


LS 


He+ 


Ballance et al. ( 


'2003) 


5 


LS 


Li+ 2 


Ballance et al. ( 


'2003) 


5 


LS 


Be+ 3 


Ballance et al. ( 


'2003) 


5 


LS 


B+4 


Ballance et al. ( 


'2003) 


5 


LS 


C +5 


Ballance et al. ( 


'2003) 


5 


LS 


o+ 7 


Ballance et al. I 


'2003) 


5 


LS 


Ne+ 9 


Ballance et al. ( 


'2003) 


6 


LS 


Fe +25 


Ballance et al. ( 


'2003) 


5 


IC 



tions of IS anrpson et all (|1983l ) for the remainder. The 
IKisielius et al.l (|1996f) data were only calculated for He 
and Fe; values for intermediate ions were obtained by 
scaling with Z . The maxmimum n-shell included in these 
data sets was n=5. 

We have incorporated a new set of i?-matrix results, 
obtained using the RMPS approach. When sophisticated 
techniques are applied to the H-like systems the results 
are not always published, as they are relatively straight- 
forward calculations used as tests of the technique itself. 
As a result, i?-matrix results could only be obtained for 
a few ions of astrophysical relevance. The ions for which 
i?-matrix data were available and their sources are listed 
in Table [2] 

For most of these data sets, calculations have been con- 
ducted in LS coupling. The resolution of these calcu- 
lations stops at the term (e.g. 2 P), and not the level 
which includes fine structure splitting by total orbital 
angular momentum J (e.g. 2 Pi, 2 Pz). Due to the small 
J-splitting of the energy levels, term resolution is often 
adequate for hydrogen-like ions. For example, the split- 
ting of the n = 2 energy levels of Ne X AE re 0.5eV 
at E re IkeV. This would therefore require a spectral 
resolution of E/AE > 2000 to resolve, while the Chan- 
dra HETG has only E/AE = 1000. The resolution 
required to observe the lines varies as approximately 
E/AE oc z -2 , where z is the atomic number of the ele- 
ment. At higher z this splitting becomes apparent - for 
Fe XXVI, E/AE re 328, which is observable in HETG 
spectra. This, in addition to the fact that AtomDB is 
used not only in analyzing existing mission data, but also 
for planning new missions which may require higher res- 
olution spectra, requires all the atomic data in AtomDB 
to be resolved into its fine structure levels if possible. 

We have therefore adopted a statistical splitting 
method for the hydrogenic ions to apply the fine- 
structure splitting: collisional excitation from ns (e.g. Is, 
2s, 3s...) levels is taken from the term resolved i?-matrix 
results: the rates are then split by statistical weights be- 
tween J-resolved levels. This method is only possible for 
the ns levels, yet this is largely sufficient for our needs 
since the ground and first excited state are ns levels, and 
these levels are the dominant sources for excitation to 
higher levels and therefore the resulting spectra. This 
method was also adopted by the CHIANTI database. 

For each H-like ion, the data available in the literature 
was expanded upon using a combination of FAC and Au- 
TOSTRUCTURE. Energy levels and transition rates were 
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calculated for and between each level from n = 1, I = 
up to n — 10, Z = 9, for a total of 99 levels, using 
Autostructure. Energy levels calculated in this way 
agreed with values from NIST to within 0.05%. Since ob- 
served wavelengths are stored separately from theoretical 
wavelengths in AtomDB and are used for spectral calcu- 
lations in preference to theoretical values where possible, 
these were not corrected further during these calculation. 
Transition probabilities were calculated for all transitions 
up to the electric octupole and magnetic quadrupole. 

Distorted wave Maxwell-averaged electron-impact col- 
lision strengths between each of these levels were cal- 
culated using FAC with the energy levels calculated by 
Autostructure. These "fill in" data sets were used in 
three ways: 

1. To provide the collisional excitation and radiative 
rates for transitions where 10 > n upper > 5 (for 



Ne +9 ,10 > n-u 



> 6) 



2. To provide the collisional excitation and radiative 
rates between nl levels, where li ower > 

3. For those ions not included in Tabled to provide 
the full range of excitation and radiative rates for 
all transitions. 

The nl level transitions are weak and are not major 
contributors to the overall spectra, however they are in- 
cluded for completeness. Also for completeness, for all 
ions not mentioned in Table [2] between H and Kr, the 
data obtained from this method were included in the 
database. 

3.4. He-like ions 

For the He-like ions, data have been incorporated from 
Intermediate Coupling Frame Transformation (ICFT) R- 
matrix calculations provided by A. Whiteford (2005, pri- 
vate communication) for all elements from C 5+ to Kr 35+ . 
In this work, structure data is again generated by Au- 
tostructure, with the ICFT method being used for 
the collision strengths. A complete description of these 
calculations has not been published in the literature, but 
the data is available at the website of the UK APAP col- 
la boration^!, and the tech nique is outlined by the author 
in lWhitefordlFall (l200l . 

Since this data is an unpublished collection, we have 
compared it with the fully relativistic .B-m atrix calcu- 
lation of collisional excitation for Ne ix by (jChen et al.1 
2006). The Maxwell -averaged collision strengths were 
found to agree with each other to within 5% for all tran- 
sitions, with the exception of those from the Is 2 1 So — > 
Is 1 ns 1 1 So, where the agreement is still within 10%. A 
best case and worst case scenario comparison is shown 
in Figure [TJ Al so included in this figu re are the collision 
strengths from ISampson et all (|1983l ). which were used 
in AtomDB vl.3.1. These data were obtained from non- 
relativistic Coulomb Born Exchange calculations, and 
only at 9 energies such that resonance effects are not 
included. 

There are two widely used diagnostic line ratios which 
arise from the n = 2 — >• 1 transitions of He-like sys- 
tems: the electron temperature sensitive G-ratio, and 

2 http : //amdpp . phys . strath . ac . uk/UK_APAP/DATA/adf 04/hel ike 
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Fig. 1. — The Maxwell-averaged collision strengths for two tran- 
sitions of Ne IX. Left Is 2 1 So — > Is 1 2s 1 1 So, and right Is 2 1 Sp — » 
Is 1 2d 1 1 P 1 . The dotted line (SGC) is the data of Sampson et al. 
l(T983h . used in AtomDB vl.3.1, the solid line (ADW) denotes the 
data of A. Whiteford (2003, private communication), the dash-dot 
line (GXC) is that of Chen et al. (2006|). The lower panel shows 
the G-ratio for both all three calculations, including recombina- 
tion cascad e s. Th e squares are the values quoted for G-given in 
IChen et al.l H2006T ), our recalculation differs slightly due to differ- 
ent handling of the recombination and cascades. 



the e lectron density sensitive R-ratio dGabriel fc Jordan! 
11969ft . In the notation of iGabriel fc Jordan! (|1969ftT 
these lines of interest are: the resonance line, w, 
(Is 2p x Pi — > ground); the intercombination lines, x and 
y, (Is 2p 3 P[i.2] - > ground); and the forbidden line, z, 
(Is 2s 3 Si — > ground). The G-ratio, (x + y + z)/w de- 
creases as the electron temperature, T ei increases, allow- 
ing a this line ratios to be used to measure T e . The de- 
crease is driven by a number of factors: the excitation of 
the forbidden level has strong low energy resonances, and 
is therefore more pronounced at low temperature. In ad- 
dition, the most prominent collisional excitation process 
changes: to excite the triplet levels of the x, y and z lines, 
the incident electron must exchange with an electron in 
the target ion, while the singlet level can be populated 
directly by excitation. The cross section for the former 
process is larger at lower energy and declines at higher 
energies while the reverse is true for direct excitation. 
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The resonance line does not completely dominate due to 
recombination and direct excitation into higher energy 
levels, and the subsequent cascade to the n = 2 shell, 
which populates them broadly in line with their statis- 
tical weights. This in turn prevents the G-ratio from 
falling significantly below 0.5. 

The R-ratio, z/(x + y), is sensitive to the electron den- 
sity, N e . At very low densities, electrons entering the 
metastable Is 2s 3 Si level decay to the ground state de- 
spite the state's very long lifetime, as there is no alterna- 
tive lower energy state to decay to. At higher densities, 
electron collisions depopulate this state, leading to a rela- 
tive decline of the forbidden line. This makes the R-ratio 
density sensitive, while also being relatively insensitive to 
temperature. 

In the lower panel of Figure [T] we compare the G-ratio 
from the different Ne IX calculations for electron impact 
excitation. The two i?-matrix methods are very similar 
in their results, showing that relativistic corrections are 
not large for Ne ix. The similarity between the fully- and 
partially-relativistic i?-matrix data is not unexpected for 
the Ne-like ions. When approaching heavier, more highly 
ionized ions (e.g. Fe xxv), relativistic effects will become 
m ore significan t and a larger discrepancy is expected. 
In lChenl (|2008f >. the authors compare the ir results with 
EBIT measurements from lWargelinl (1993) and find very 
good ag reement. This d iffers significantly from the re- 
sults of I Sampson et al.1 (|1983| ) . We have incorporated 
the data of A. Whiteford for the He-like isoelectronic se- 
quence into AtomDB v2.0.2: it provides a high-quality 
uniform data set covering a wide range of ions. 

As with the H-like ions we have extended the range of 
these calculations using FAC to obtain collision strengths 
for collisions with 5 < n < 10. This was done using 
exactly the same method as for the H-like ions. 

3.5. Iron L-shell Data 

The data for the iron L-shell ions (Fe +16 to Fe +23 ) 
in At omDB vl.3.1 is taken from a series of HUL- 
LAC (|Bar-Shalom et al.ll200l calculations by D. Liedahl 
(1997, private communication). While these data 
have never been explicitly published, the data for 
Fe + W to Fe +23 are descri bed in iLiedahl et all (|1995f ) 
and iWargelin et al.l (fl998). In summary, the struc- 
ture was calculated using the relativis tic, muticonfig- 
uration parametric potential method (Klapisch 1971; 
iKlapisch et al.l Il977t ) , while the electron collision data 
are semi-relativistic distorted wave calculations, and in- 
clude excitation and radiative decay from up to the n = 5 
to n = 7 shell, depending on the ion. 

Significant advances have been made in collisional ex- 
citation calculations since the release of AtomDB vl.3.1. 
These ad vances are large l y due to the Iron Project, ini- 
tiated bv iHummer et ail (|1993f l. We have taken advan- 
tage of this data to upgrade the collision strength data 
for these ions to i?-matrix calculations. Table [3] lists 
the references for each ion in the Fe L-shell series. These 
new calculations include, in general, fewer configurations 
than the existing HULLAC data; as a result the newer 
data have been merged with the existing data, updat- 
ing or adding as appropriate. Much of the older data 
remains for transitions not included in the newer data. 
Table |4] shows the configurations included for each ion in 
both the new and old datasets; unless otherwise noted, a 



TABLE 3 

The data sources used for iron L-shell 
ions, and the maximum n shell included 

FOR EACH. 



Ion Reference n mas 

Fe+ 16 Loch et al. (2006') 4 

Fe+ 17 Witthoeft et al. (2007b) 4 

Fe+ 18 Butler fe Badnell (20081 4 

Fe+ 19 Witthoeft et al. (2007a) 4 

Fe+ 20 Badnell fe Griffin (2001) 4 

Fe+ 21 Badnell et al. (2001) 4 

Fe+ 22 Chidichimo et al. (2005 ) a 4 

Fe+ 23 Whiteford et al. (2002) 3 

a as amended bv lDel Zanna et al.l (2005) 



level included in the new i?-matrix data is also included 
in the old HULLAC data. In general the HULLAC data 
includes n shells two or three higher than those in the 
i?-matrix data. 

Merging the two data sets requires a method for match- 
ing levels between the two different data sets. This is 
complicated by the different coupling schemes reported 
for both data sets: the Liedahl HULLAC data is jj cou- 
pled, while the new i?-matrix data is LS coupled. Levels 
have been matched by grouping levels with the same elec- 
tron shell occupancies and J quantum number, sorting 
these into energy order, and linking it to its counterpart 
in the other data set after the performing the same pro- 
cess. Where level energies disagree, the i?-matrix calcu- 
lations have been used. Spot checks on many strong lines 
show that these have been correctly matched, although it 
is difficult to check for the vast number of levels without 
strong emission. 

Due to its abundance in the Universe and its strong 
emission over a wide range of wavelengths, line ratio diag- 
nostics based on Fe xvn should be very powerful. There 
are, however, significant differences in the values of these 
line ratios in the literature, making their use problematic. 

The most discussed diagno stic line ratios, following no- 
tation from iParkinsonl (|1973f ). are the 3C/3D ratios and 
the 3s/3C ratios, where 3C = 2p 5 3d 1 1 Pi — > ground, 
15.014A; 3D = 2p 5 3d 1 3 Di -> ground, 15.26lA; and 3s 
is the sum of the 2p 5 3s 1 — > ground transitions. The 
3C/3D line ratio has long been proposed as a measure of 
the opacity of the plasma, since the 3C line is sensitive 
to the opacity of the plasma and the 3D line is less so: 
comparing the observed ratio with the prediction for an 
optically thin plasma can be used to obtain the optical 
depth of the plasma. 

Several groups have measured these line ratios 
using EBIT e xperim ents, includi ng those at NIST 
(ILaming et all I2000D a nd LLNL (jBrown et al.l I1998L 
120011: iBeiersdorfer et al.ll2002f) . These results have been 
inconsistent with each other, with cal culated value s 
of the rat i os for thes e lines, such as iChenl {2008), 
ILoch et alJ (|2006l ) and ILaming et aT] (12000D. and with 
astrophysically observed values (|Brickhouse fe Schmeld 
I2006T ) . Agreement was frequently no better than a factor 
of 2 between the differ ent methods. 

Recently the work of IGillaspv et al.l (|2011f ) appears to 
have shown that with careful consideration of the EBIT 
plasma and instrumental effects such as window absorp- 
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TABLE 4 

The configurations included in the iron L-shell datasets listed in Table [3] and the 
additional configurations included in the atomdb vl.3.1 hullac datasets. 



Ion 


Configurations (New data) 


/iQaitionai oonngurations (^riuijijrt-^ aaxaj 


Fe +16 


ls 2 2s 2 2p 6 , ls 2 2s 2 2p 5 31 1 , ls 2 2s 1 2p 6 31 1 , 


ls 2 2s 2 2p 5 61 1 , ls 2 2s 2 2p 5 71 1 , ls 2 2s 1 2p 6 51 1 , 




ls 2 2s 2 2p 5 41 1 , ls 2 2s 2 2p 5 51 1 , ls 2 2s 1 2p 6 41 1 


ls 2 2s 1 2p 6 61 1 , ls 2 2s 1 2p 6 71 1 


Fe+ 17 


ls 2 2s 2 2p 5 , ls 2 2s 1 2p 6 , ls 2 2s 2 2p 4 31 1 , 


ls 2 2p 6 31 1 , ls 2 2s 1 2p 5 41 1 , ls 2 2s 2 2p 4 51 1 , 




ls 2 2s 1 2p 5 31 1 , ls 2 2s 2 2p 4 41 1 


ls 2 2s 1 2p 5 51 1 , ls 2 2p 6 41 1 , ls 2 2p 6 51 1 


Fe+W 


ls 2 2s 2 2p 4 , ls 2 2s 1 2p 5 , ls 2 2p 6 , 


ls 2 2s 1 2p 4 41 1 , ls 2 2s 2 2p 3 51 1 , ls 2 2p 5 41 1 , 




ls 2 2s 2 2p 3 31 1 , ls 2 2s 1 2p 4 31 1 , ls 2 2p 5 31 1 , 


ls 2 2s 1 2p 4 51 1 , ls 2 2p 5 51 1 




ls 2 2s 2 2p 3 41 1 




Fe+W 


ls 2 2s 2 2p 3 , ls 2 2s 1 2p 4 , ls 2 2p 5 , 


ls 2 2p 4 31 1 , ls 2 2s 1 2p 3 41 1 , ls 2 2s 2 2p 2 51 1 , 




ls 2 2s 2 2p 2 31 1 , ls 2 2s 1 2p 3 31 1 , ls 2 2s 2 2p 2 41 1 


ls 2 2p 4 4l\ ls 2 2s 1 2p 3 51 1 


Fe+ 2 


ls 2 2s 2 2p 2 , ls 2 2s 1 2p 3 , ls 2 2p 4 , 


ls 2 2p 3 3p 1 , ls 2 2p 3 3d 1 , ls 2 2s 1 2p 2 41 1 , 




ls 2 2s 2 2p 1 31 1 , ls 2 2s 1 2p 2 31 1 , ls 2 2p 3 3s 1 , 


ls 2 2s 2 2p 1 51 1 , ls 2 2p 3 41 1 , ls 2 2s 1 2p 2 51 1 




ls 2 2s 2 2p 1 41 1 




Fc+ 21 


ls 2 2s 2 2p 1 , ls 2 2s 1 2p 2 , ls 2 2p 3 , 


ls 2 2s 2 51 1 , ls 2 2s 1 2p 1 51 1 , ls 2 2s 2 61 1 , 




ls 2 2s 2 31 1 , ls 2 2s 1 2p 1 31 1 , ls 2 2s 2 41 1 , 


ls 2 2s 1 2p 1 61 1 




ls 2 2s 1 2p 1 41 1 , ls 2 2p 2 31 la 




Fe+22 


ls 2 2s 2 , ls 2 2s 1 2p 1 , ls 2 2p 2 , 


ls 2 2s 1 51 1 , ls 2 2p 1 51 1 , ls 2 2s 1 61 1 , 




ls 2 2s 1 31 1 , ls 2 2p 1 31 1 , ls 2 2s 1 41 1 , 


ls 2 2p 1 61 1 , ls 2 2s 1 71 1 , ls 2 2p 1 71 1 




ls 2 2p 1 41 1 




Fe+23 


ls 2 21 1 , ls^l 1 , ls 2 41 1 , 


ls 2 61 1 , Is 2 ?! 1 




ls 2 51 1 , ls 1 2s 2a ls 1 2s 1 2p lB , 






ls 1 2p 2a , ls 1 2s 1 31 la , ls 1 2p 1 31 la 





a These configurations were included in the new datasets but were not in the HULLAC 
dataset from AtomDB vl.3.1 



tion, as well as revisiting RR recombination and its pop- 
ulating effect on upper levels, that observed values can 
agree to within 5 to 20% of t he calculated values. We 
show in Figured adapted from[Chei (|2008l ). the 3C/3D 
line ratio in iron from a variety of astrophysical objects. 
Also shown are the 3C/3D ratios form the Raymond- 
Smith co de, from Ato mDB vl.3.1, from AtomDB v2.0.2, 
and from[Chen| (12001 . 

In this c ase, the AtomDB v2.0.2 data is t hat of 
ILoch et ail (|2006j ) . We note that the values from IChenl 
(2008) are closer to the EBIT measurements and to the 
mean of the observed values, and therefore that their 
line ratio values are more likely accurate. However, 
this calculation is incomplete for the remaining strong 
transitions of Fe XVII, and as individual line ratios are 
a complex result of many different excitation and de- 
excitation processes, this data cann ot be included in 
AtomDB. We have therefore used the ILoch et~aH (|2006f) 
data in AtomDB v2.0.2, providing a significant improve- 
ment over the data from vl.3.1, but this will be revisited 
when further data becomes available; for now an uncer- 
tainty of up to 2 0% should be assumed w hen using values. 

As noted by IDoron fe Beharl (|2002f ) . the addition of 
the state-selective recombination rates leads to some sup- 
pression of the 3C/3D line ratio at high temperature. 
They approached this issue by constructing a three ion 
model to include ionization and recombination, and our 
results here do not conflict with their findings. 

3.6. Other Ions 

For many other ions, especially those that are not no- 
tably strong emitters in the X-ra y region, AtomDB vl.3.1 
used data from CHIANTI v2.0 (|Landi et al.|[T999h . We 
have updated the wavelengths, transition probabilities 
and collision strengths of all data for all ions excluding 
the H-like & He-like sequences, and the Ni and Fe L- 
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Fig. 2. — The Fe XVII 3C/3D line ratio from various calcula- 
tions, compared with observed data. Solid line: AtomDB v2.0.2; 
dash-dotted line: AtomDB v2.0.2 excluding state-selective recom- 
bination; dotted line : AtomDB vl.3.1; small dashed line: th e Ray- 
mond Smith model (Raymond & Smith 1977); large dashes: IChenl 
(2008). Various observed values are also taken from a variety of 
Chandra and XMM-Newton obser vations (solid and open circles, 
respectively). Figure adapted from Chen (2008). 



shell ions to now incorporate the CHIANTI v6.0.1 data 
(pere et al.l 12009ft . which was the current version at the 
time of preparation of this data set. The Ni L-shell data 
remains that of the previous version of AtomDB, which 
is scaled from the Fe data of Liedahl. We will aim to 
address this in a future update. 



3.7. Minor Changes 

In addition to these major updates to the database, 
several minor corrections have been made which have lit- 
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tie or no effect on the resulting spectra from AtomDB but 
which do affect the database. These include the correct 
use of L and S quantum numbers for all levels where LS 
coupling is used; consistent formatting of configuration 
strings, minor corrections to the H-like two-photon tran- 
sition rate from the 2s 1 — > Is 1 level, including adding the 
rate for hydrogen and reducing it for all other H-like ions. 
None of these changes have a noticeable effect on the 
spectra from AtomDB v2.0.2 compared with AtomDB 
vl.3.1. 

4. RESULTS 

4.1. Effects of Improved Ionization and Recombination 

Rates 

Changes in the ionization balance are most significant 
for heavier eleme nts. Detailed information c an be found 
in the papers of iBrvans et all (pOOl [2009h . Figure 02 
shows the change in the ionization balance for iron, in- 
cluding the recombination rates for Fe XIX as specified 
in Section 13.11 Significant changes can be seen around 
Fe xvn, in some cases leading to changes in fractional 
abundance of up to 30%. 

The most readily apparent spectral changes between 
AtomDB vl.3.1 and v2.0.2 are due to the change in the 
ionization balance data. The shift in the ionization bal- 
ance leads to stronger emission from many of the iron 
L-shell ions at temperatures between 10 6 and 10 7 K. In 
the lower part of Figure 02 the line emissivities from 
AtomDB vl.3.1 are compared to that of AtomDB v2.0.2 
for a 1 keV optically-thin, collisionally-ionized plasma, 
with N e = lcm" 3 . Also, a combinati on of the v2.0.2 
data w ith the vl.3.1 ionization balance oflMazzott a et all 
(1998) is shown, to highlight the effects of the new ion- 
ization balance as opposed to all the other new data. As 
can be seen, the majority of the change in emission is 
due to the ionization balance, which leads to increased 
Fe xvni and Fe xix emission in the 0.7 to 1 keV band. 
The continuum emission in all three cases changes very 
little, and therefore is not plotted in detail. 

4.2. Cooling Function 

Figure 0] shows the total power radiated from all ele- 
ments given a solar photo spheric elemental abundance 
(|Anders fe Grevesse] I1989D and an electron density of 
1 cm~ 3 . The top figure shows a comparison between 
the total emission from the old and new versions of 
the database: as can be seen, there are significant (> 
10%) increases in radiation in three different tempera- 
ture ranges. These are due to changes in the ionization 
balance for oxygen, neon and iron respectively: in the 
new data a higher fraction of less-ionized ions is present. 
These ions radiate more strongly in higher temperature 
regions, in part due to increased recombination-cascade 
emission. Figure |4] also shows the radiation separated by 
element as a function of temperature. 

4.3. Higher n-shell Effects 

The inclusion of higher n-levels for H- and He-like sys- 
tems has two effects: in addition to new lines formed 
from these levels, cascade effects to lower n levels can 
also alter line intensities. Figure 03 shows the emissivi- 
ties of lines of Ne X with the higher n-shell lines appar- 
ent and demonstrates the possibility of either observing 
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Fig. 3. — Top: T he ionization ba l ance f or iron above 10 6 K. The 
dashed line is the Mazzotta ct al. (1998) data use d in AtomDB 
vl.3.1 , the s olid line is the newer data mostly based on lBrvans et al.l 
(2006, 2009) included in v2.0.2. Significant chan ses are found in 
the 10 6 — 10 7 K range. Bottom: A comparison of the continuum 
and line emission from vl.3.1 and v2.0.2 of AtomDB, in a IkeV 
optically thin plasma with solar abundances |Andcrs & Grevesse] 
1989). Lines have been broadened by convolution with Gaussians, 
with FWHM=0.0235 keV. Also shown is the line emission using 
all of the atomic data from AtomDB v2.0.2, but the ionization 
balance from vl.3.1. The ionization balance is the dominant cause 
of changes in the spectrum, with more Fe+ 17 and Fe+ 18 present in 
version. 



these lines, or, more likely, having these lines blend with 
other lines from Mg XI. 

Figure [6] shows the Mg xi bandpass (9-10A) of 
the Chandra HETG spectrum from TW Hydrae 
(|Brickhouse et al.ll2010h . compared to a simple optically- 
thin, collisionally ionized plasma spectrum created us- 
ing AtomDB v2.0.2. Line features are broadened by 
convolution with Gaussians with FWHM=0.010A. Con- 
tribution from Ne X to the spectrum is highlighte d in 
red. Solar abundances of I Anders fe Grevesse! (|1989h are 
used throughout, except for Ne (=1.23 solar) and Mg 
(=0.18 solar). These va lues are taken from Model A in 
iBrickhouse et all ([2010), as are the plasma conditions of 



5.8 x 10 12 cm 



The inclusion 



T e = 3.58 x 10 6 K, 
of the n > 5 levels leads to several extra lines being ap- 
parent in the spectrum, potentially affecting line ratio 
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Fig. 5. — Top: the line emissivities for Ne X at 10 7 K, showing 
the old data (dashed line) and the new (solid). The lines from the 
high-n cascade can be seen clearly. Bottom: A simulated spectrum 
of Ne X from a lOOks Chandra HETG observation, also at 10 7 K. 



Fig. 4. — Top: the cooling function in t he 0.001 to lOOkeV 
band from all ions, assuming solar abundances ( Anders & Grevessel 
1989), comparing the old (vl.3.1, dashed) and new (v2.0.2, solid) 
AtomDB. Bottom The cooling function for AtomDB v2.0.2, broken 
down by element. 



estimates depending on the resolution of the instrument 
in use. 

4.4. Diagnostic Line Ratios 
4.4.1. He-like ratios 

The effect on the G-ratios (Section l3.4l) of the new data 
used in AtomDB v2.0.2 can be seen in Figure[7J For most 
ions the G-ratio substantially increases, in some cases 
leading to a doubling of the implied temperature. At the 
peak temperatures for each ion, this effect is attributed 
largely to an increase in the collisional excitation rate to 
the forbidden and intercombination transition upper lev- 
els: this has led to, for example, a 30% increase in both 
the forbidden and intercombination line emission for Si 
XIII. At slightly higher temperatures, above which the 
n = 2 — > 1 complex emission has peaked, the recombi- 
nation and cascade processes have the most significant 
effect. 

Comparisons of the temperature implied by the G-ratio 
with that implied by the ratio of the Ly-a/w, where w is 
the He-like resonance line are sometimes used to investi- 
gate plasma ionization states. The two line ratios mea- 
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Wavelength [A] 

Fig . 6.— Top: Summed C handra HETG spectrum of TW Hy- 
drae (Brickhousc ct al. 2010) in the Mg XI bandpass. Bottom: A 
simulated AtomDB spectrum with Gaussian broadening applied to 
each line. Contributions from Ne X are shaded in red. 



sure slightly different temperatures. The G-ratio mea- 
sures directly the electron temperature of the plasma, 
due to the fast electron collision rates. The Ly-a/w ra- 
tio, however, primarily depends on the ionization bal- 
ance, and therefore can be used to indicate the ioniza- 
tion temperature of the plasma. If these two diagnos- 
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Fig. 7. — The temperature sensitive G-ratio for selected He-like 
ions, for both the old (vl.3.1, dashed) and new (v2.0.2, solid) 
AtomDB. 



tics disagree about the electron temperature, then the 
plasma is either over- or under-ionized; this is possible 
since ionization and recombination occur on much slower 
time-scales than collisional excitation. For a single tem- 
perature plasma in ionization equilibrium they will be 
equal; however, in an object with a broad thermal dis- 
tribution, the hy-a/w ratio can be elevated due to Ly-a 
e mission from warm er plasma. 

iNess et all (|2003t ) encountered this problem when com- 
paring the G-ratio for Ne IX with the emission measure 
obtained from iron L-shell ions in Capella: the emission 
measure peaks at aro und 6MK, but the G-ratio implied 
a T e of around 2MK. iTesta et all (|2004f) noted that in 
Chandra observations of a variety of stars, the temper- 
ature was consistently underestimated by the G-ratios, 
sometimes by up to a factor of four. In Figure [8] we 
repeat this comparison, using the data from AtomDB 
v2.0.2, and the line fluxes from Tables 6 and 7 of their 
paper. Using the AtomDB vl.3.1 data, we have been 
unable to reproduce the exact temperatures from their 
p aper, but the same tren d is clearly evident . 

ITesta et al.l (|2004f ) and INess et all (|2003l) discuss the 
causes of some of these discrepancies. The differential 
emission measure (DEM) of the stars in the survey is gen- 
erally peaked around 10 7 K (see references in ITesta et al.l 
2004) . The temperatures of peak emissivity of O VII and 
Mg xii are » 2 x 10 6 K and « 6.4 x 10 6 K respectively. By 
10 7 K, the emissivity of the O VII has dropped two orders 
of magnitude from its peak. Therefore we would expect 
the G-ratio from oxygen to be slightly higher than the 
peak emissivity at 2 x 10 6 K, but not as high as 10 7 K due 
to the lack of O vn emission at this temperature. The 
Mg xi G-ratio should give a temperature closer to the 
peak of the DEM at 10 7 K as the emissivity of this ion is 
still high. Instead, they found that the G-ratio in both 
cases gave temperatures which were below the tempera- 
tures of peak emissivity, which suggested a problem in 
the atomic data. 

With the AtomDB version 2.0 data, we find that the 
G-ratio now indicates temperatures which lie on aver- 
age at or above the peak emissivity temperature for each 
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Fig. 8. — A comparison of the temperatures implied by the mea- 
sured G-ratios a nd Ly-g/resonance lines for Chandra MEG obser- 
vation data from Testa et al. (2004), using AtomDB vl.3.1 data (as 
in the original paper) and the new AtomDB v2. 0,2 data. Er r or bar s 
are based on the error in the line fluxes from ITesta et al] 1120041 ); 
X-axis errors are very small due to the very strong temperature 
dependence of the Ly-a/w ratio. 



ion (some spread remains due to the different DEMs of 
different objects). In addition, the discrepancy between 
the two temperature measurement methods has been re- 
d uced from 60% to a round 35%. 

ISmith et al.l (|2009t ) also investigated the G- and R- 
ratios for the case of Ne IX. They used AtomDB vl.3.1 
data, however they updated it t o inclu de the more recent 
collisional data of IChen et al.l (|2006D (discusse d in sec- 
tion | 3 .4I) and the recombination data from the IBadnelll 
( 2006a) collection - the same data included in the release 
of AtomDB v2.0.2. The authors found a significant vari- 
ation with temperature in the R-ratio, which is the ratio 
of the forbidden to the inter-combination lines of He-like 
ions. This temperature dependence is unusual: the G- 
ratio is typically a widely used density diagnostic with 
little temperature variation. We do not find any tem- 
perature dependence here, despite using almost identical 
data (the small differences in the collisional excitation 
calculations ha ve almost no effect ). This was traced to 
an error m the ISmith et~aT1 {2009) paper when connect- 
ing the energy levels in the recombination data to the 
excitation data. Figure UH shows the corrected R-ratio as 
a function of temperature. 

5. SUMMARY AND FUTURE PLANS 

We have presented the latest version of AtomDB, the 
first major update since 2001. Nearly every single piece 
of data in the database has been updated, with many 
new ions added. This is the result of a comprehensive 
evaluation of the previous data, and assessment of its re- 
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Fig. 9. — The density diagnostic R-ratio from AtomDB vl.3.1 
( dash-dot), t he At omDB v2.0.2 data (solid), and the data from 
ISmith ct al. (2009]) both with (dotted) and without (dashed) re- 
combination. 



placement, addressing many of the known issues in the 



previous version. New recombination data significantly 
alter the ionization balance, and therefore the emissivi- 
ties, of many lines; new collisional data for H- and He-like 
ions make significant improvements to common temper- 
ature diagnostics. 

The new data are now available on-line at 
www.atomdb.org, and are also avai lable throug h 
spectr al fitting packages such as XSPEC (jArnau 311990) . 
ISIS (|Houck fc Denicolal |2000!) and Sherpa (|Freeman| 
[200lT) FI 

There are several improvements already planned for 
the next release of AtomDB. We will include final state 
resolved DR rates for the remaining ions in the database 
for which the data exist. We will largely target non- 
equilibrium ionization plasmas, with new inner-shell ex- 
citation data, as well as fluorescence line data. Work is 
also ongoing to document and release the APEC colli- 
sional ionization code, which will allow users to generate 
the non-equilibrium higher density plasma models. 

We wish to thank Paola Testa for discussion in- 
volving the line diagnostic issues. The authors grate- 
fully acknowledge funding from NASA ADP grant 
#NNX09AC71G. 
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